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g ; ABSTRACT 

£N1 ' Microlensing is one of the promising techniques that can be used to search for extra-solar 

systems. Planet detection via microlensing is possible because the event caused by a lens system 
having a planet can produce noticeable anomalies in the lensing light curve if the source star 
passes close to the deviation region induced by the planet. Gaudi, Naber & Sackett pointed 
out that if an event is caused by a lens system containing more than two planets, all planets 
will affect the central region of the magnification pattern, and thus the existence of the multiple 
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planets can be inferred by detecting additionally deformed anomalies from intensive monitoring of 
, high magnification events. Unfortunately, this method has important limitations in identifying 

the existence of multiple planets and determining their parameters (the mass ratio and the 
instantaneous projected separation) due to the degeneracy of the resulting light curve anomalies 
from those induced by a single planet and the complexity of multiple planet lensing models. 
In this paper, we propose a new channel to search for multiple planets via microlensing. The 
method is based on the finding of Han et al. that the lensing light curve anomalies induced by 
multiple planets are well approximated by the superposition of those of the single planet systems 
where the individual planet-primary pairs act as independent lens systems. Then, if the source 
trajectory passes both of the outer deviation regions induced by the individual planets, one can 
unambiguously identify the existence of the multiple planets. We illustrate that the probability 
of successively detecting light curve anomalies induced by two Jovian-mass planets located in the 
lensing zone through this channel will be substantial. Since the individual anomalies can be well 
modeled by much simpler single planet lensing models, the proposed method has an important 
advantage of allowing one to accurately determine the parameters of the individual planets. 

Subject headings: gravitational lensing - planets and satellites: general 

1. Introduction possible because the lensing event caused by a lens 

system containing a planet can produce notice- 
Microlensing is one of the promising techniques able anomalies in the resulting light curve when 

that can be used to search for extra-solar plan- the source passes c i ose to the lens caustics, which 

ets, especially located at large distances (Perry- represents the source positions on which the lens- 

man 2000). Planet detection via microlensing is ing magnification of a point source event becomes 
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infinity (Mao & Paczyriski 1991; Gould & Loeb 
1992). 1 For a lens system with a planet, there exist 
two or three disconnected sets of caustics. Among 
them, one is located close to the primary lens (cen- 
tral caustic) and the other (s) is (are) located away 
from the primary lens (planetary caustic[s]). Ac- 
cordingly, there exist two types of planet-induced 
anomalies: one affected by the planetary caustic 
(type I anomaly) and the other affected by the cen- 
tral caustic (type II anomaly) (Han & Kim 2001). 
Compared to the frequency of type I anomalies, 
type II anomaly occurs with a relatively low fre- 
quency due to the smaller size of the central caus- 
tic compared to the corresponding planetary caus- 
tic. However, the efficiency of detecting type II 
anomalies can be high because intensive moni- 
toring is possible due to the predictable time of 
anomalies, i.e. near the peak of magnification, and 
the known type of candidate events for intensive 
follow-up monitoring, i.e. very high magnification 
events (Griest & Safizadeh 1998). 

Keeping the high efficiency of type II anomaly 
detections in mind, Gaudi, Naber & Sackett 
(1998) pointed out that if an event is caused by 
a lens system having multiple planets located in 
the lensing zone, within which the chance for the 
occurrence of planet-induced anomalies is maxi- 
mized 2 , all planets will affect the central region 
of magnification pattern, and thus the existence 
of multiple planets can be inferred by detecting 
additionally deformed anomalies from intensive 
monitoring of high magnification events. This 
method, however, has important limitations in 
identifying the existence of multiple planets and 
determining their parameters (the mass ratio and 
the instantaneous separation between each planet 
and the host star). This is because the anoma- 
lies induced by multiple planets are in many cases 
qualitatively degenerate from that induced by a 
single planet and even if the existence of the mul- 
tiple planets is known, accurate determination of 
the individual planet parameters will be difficult 
due to the complexity of multiple planet lensing 
models. 

In this paper, we propose a new channel to 
search for extra-solar systems composed of multi- 
ple planets via microlensing. The method is based 



1 See § 2 for more details about the caustics. 

2 See also § 2 for more details about the lensing zone. 



on the finding of Han et al. (2001) that the lensing 
light curve anomalies induced by multiple plan- 
ets are well approximated by the superposition 
of those of the single planet systems where the 
individual planet-primary pairs act as indepen- 
dent lens systems. Then, if the source trajectories 
passes both of the outer deviation regions around 
the planetary caustics of the individual planets, 
one can identify the existence of the multiple plan- 
ets. We illustrate that the probability of succes- 
sively detecting lensing light curve anomalies in- 
duced by two Jovian-mass planets located in the 
lensing zone through this channel will be substan- 
tial. We discuss the advantages of the proposed 
method over the previous method of monitoring 
high magnification events. 

2. Basics of Multiple Planet Lensing 

If a source located at r s on the projected plane 
of the sky is lensed by a lens system composed of 
iV-point masses, where the individual components' 
masses and locations are m, and 1$, the positions 
of the resulting images, r, are obtained by solving 
the lens equation of the form 



r s — r - 



where m = rrii is the total mass of the lens 
system and 6*e is the angular Einstein ring radius. 
The angular Einstein ring radius is related to the 
physical parameters of the lens system by 

flE=l /^(^-^V /2 , (2) 



N 

q2 \ - rrii 
m 



i=l 



k\ 2 ' 



(1) 



D„ 



where D a i and D os are the distances to the lens 
and the source from the observer, respectively. 
Since the lensing process conserves the surface 
brightness of the source, the magnification of each 
image is simply given by the surface area ratio 
between the image and the unamplified source, 
and mathematically its value corresponds to the 
inverse of the Jacobian of the lens equation evalu- 
ated at the image position rj ; 

dr. 
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The total magnification is then given by the sum of 
the individual images' magnifications, i.e. A tot — 



J2f* Aj, where Ni is the total number of images. 



2 



For a single point-mass lens (N = 1), there are 
two solutions for the lens equation (Nj = 2) and 
the resulting total magnification is expressed in a 
simple analytical form of 



u? + 2 



(4) 



where u s = (r s — 1)/6e is the dimensionless lens- 
source separation vector normalized by 9e- For a 
rectilinear lens-source transverse motion, the sep- 
aration vector is related to the single lensing pa- 
rameters by 



t-tp 
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where ts represents the time required for the 
source to transit 9e (Einstein time scale), [3 is the 
closest lens-source separation in units of # E (im- 
pact parameter), to is the time at that moment, 
and the unit vectors x and y are parallel with and 
normal to the direction of the relative lens-source 
transverse motion, respectively. 

The lens system with planets is described by the 
formalism of multiple lensing with very low mass 
companions. For a binary lens system (N = 2), 
the lens equation becomes a fifth order polyno- 
mial in r and the positions of the individual images 
are obtained by numerically solving the equation 
(Witt 1990). This yields three or five solutions 
depending on the source position with respect to 
the lenses. The main new feature of a multiple 
lens system is the formation of caustics. For a bi- 
nary lens, the set of caustics form several discon- 
nected closed curves. The number and locations 
of the caustic curves are dependent on the separa- 
tion between the planet and the primary lens. If 
the separation is larger than 9e, there exist two 
sets of caustics and one is the central caustic lo- 
cated near the primary lens and the other is the 
planetary caustic located away from the primary 
on the planet side with respect to the primary. 
The lens system having a planet with a separation 
smaller than 9e also has a single central caustic, 
but has two planetary caustics located on the op- 
posite side of the planet. The planetary caustic(s) 
is (are) located within the Einstein ring when the 
projected separation between the planet and the 
primary (normalized by 9e) is in the lensing zone 
of 0.618 < b < 1.618. Since the sizes of both the 



central and planetary caustics are maximized at 
around this separation, planet detection probabil- 
ity is also maximized for systems having planets 
located in this lensing zone (Gould & Loeb 1992). 

For a lens system with triple lenses (N = 3), 
the lens equation becomes a tenth order polyno- 
mial and it is still numerically solvable. For these 
lens systems, there are a maximum of ten images 
and a minimum of four images, and the number 
of images changes by a multiple of two as the 
source crosses a caustic. Unlike the caustics of a 
binary lens system forming separate sets of closed 
curves, those of a triple lens system can exhibit 
self-intersection and nesting. 

For lens systems with even larger number of 
lenses (N > 4) , numerically solving the lens equa- 
tion becomes very difficult. However, one can still 
obtain the magnification patterns of these lens sys- 
tems by using the inverse ray-shooting technique 
(Schneider & Weiss 1986; Kayser, Refsdal & Sta- 
bcll 1986; Wambsganss, Paczyhski & Schneider 
1990). In this method, a large number of light 
rays are uniformly shot backwards from the ob- 
server plane through the lens plane and then col- 
lected (binned) in the source plane. Then, the 
magnification pattern is obtained by the ratio of 
the surface brightness (i.e., the number of rays per 
unit area) on the source plane to that on the ob- 
server plane. Then the light curve resulting from 
a particular source trajectory corresponds to the 
one-dimensional cut through the constructed mag- 
nification pattern. This technique has an advan- 
tage of allowing one to construct magnification 
patterns regardless of the number of lenses, but it 
has a disadvantage of requiring large computation 
time to obtain smooth magnification patterns. 

3. A New Channel 

Recently, Han et al. (2001) showed that the 
light curve anomalies induced by multiple plan- 
ets are well approximated by the superposition 
of those of the single planet lens systems where 
the individual planet-primary pairs act as inde- 
pendent lens systems. As pointed out by Gaudi et 
al. (1998), then, the anomaly pattern in the central 
region caused by one of the planets can be signifi- 
cantly affected by the existence of other planet (s) 
because the central deviation regions caused by 
the individual planets occur in this same region. 
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Fig. 1. — The map of magnification excess of an example lens system composed of two planets with mass 
ratios of 0.003 and 0.001, which are comparable to those of Jupiter-mass and Saturn-mass planets around 
a ~ 0.3 Mq star, respectively. The coordinates (£,??) are set so the the primary is at the origin and £ axis 
is parallel with the axis connecting the primary and the Jupiter-mass planet. All lengths are scaled by the 
angular Einstein ring radius. The projected separations (normalized by #e) of the individual planets from 
the host star are bj = 1.2 and b s = 0.8, respectively. The angle between the the position vectors to the 
planets from the host star is a = 135°. Contours are drawn at the levels of e = —0.1, —0.05, 0.05, and 0.1 
and the regions of positive deviations are distinguished by grey scales. The small figures drawn by thick 
solid lines are the caustics and the dashed circle is the Einstein ring. The region enclosed by solid circle 
around the primary with a radius u s = 0.15 represents the central deviation region, which is excluded in 
our probability estimation of detecting both planets. The straight lines with arrows represent the source 
trajectories of several example events for which one can detect both (the trajectory designated by a number 
'1'), one ('3' and '4'), and none ('2') of the planets, respectively. 
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However, the anomaly pattern in the outer devia- 
tion region caused by each planet is hardly affected 
by other planet (s) because the outer deviation re- 
gions of the individual planets occur, in general, 
at different locations. 

The fact that the anomaly patterns in the outer 
deviation regions are scarcely affected by other 
planet (s) provides a new channel to search for mul- 
tiple planets and to determine their parameters. 
In this channel, multiple planets are detected when 
the source trajectory passes both of the outer devi- 
ation regions around the planetary caustics of the 
individual planets. The greatest advantage of this 
method over the method of monitoring high mag- 
nification events is that one can unambiguously 
identify the existence of multiple planets from the 
unique pattern of consecutive anomalies and accu- 
rately determine the planet parameters due to the 
applicability of much simpler single planet lensing 
models to the individual anomalies. 

4. Efficiency of the New Channel 

In this section, we illustrate that the efficiency 
of detecting multiple planets through the new 
channel will be substantial. To demonstrate this, 
we estimate the probability of successively detect- 
ing lensing light curve anomalies induced by mul- 
tiple planets. The probability is determined for a 
lens system having two Jovian-mass planets (i.e. 
Jupiter and Saturn) orbiting a host star with a 
mass of ~ 0.3 Mq. Then, the mass ratios of the 
planets (to that of the host star) are qj — 0.003 
and qs = 0.001, for the Jupiter-mass and the 
Saturn-mass planets, respectively. If the lens sys- 
tem is located located at D a i ~ 6 kpc (and with a 
source at D os <~ 8 kpc) and the component planets 
have intrinsic orbital separations similar to those 
of Jupiter and Saturn of our solar system, the or- 
bital separations (normalized by 9e) of the planets 
correspond to aj ~ 2.7 and as ~ 5.0, respectively. 

To estimate the probability, we first construct 
maps of magnification excess. The magnification 
excess is defined by 



where A and Aq represent the magnifications ex- 
pected with and without the planets, respectively. 
Figure 1 shows the constructed map when the pro- 
jected separations of the individual planets are 



bj = 1.2 and b$ = 0.8 and the angle between the 
position vectors to the planets fro the host star 
(orientation angle) is a = 135°. We note that the 
projected separation, b, is related to the intrinsic 
orbital separation, a, by 



b = ay cos 2 (f> + sin 2 (f> cos i, (7) 

where <j) is the phase angle and i is the inclination 
of the orbital plane. 

Once the map is constructed, the probability is 
estimated by computing the ratio of events whose 
source trajectories pass both of the outer devia- 
tion regions induced by the individual planets (for 
example, the event resulting from the source tra- 
jectory designated by '1' in Fig. 1) during the mea- 
surements to the total number of trial events. For 
given lens positions, the source trajectory orien- 
tations (with respect to the Jupiter-primary axis) 
and the impact parameter (with respect to the pri- 
mary) of the trial events are randomly selected in 
the ranges of < 9 < 2n and < (3 < 1.0, respec- 
tively. Measurements for each event are assumed 
to be carried out during — tE < t t, s < tE with a 
frequency of 10 times/night, which corresponds to 
that of the current microlensing followup observa- 
tions (Rhie et al. 1999; Albrow et al. 1998; Bond 
et al. 2001). We consider each planet is detected if 
excesses greater than a threshold value of e t h = 5% 
are consecutively detected more than 5 times dur- 
ing the measurements. Since we are interested in 
anomalies occurred only in the outer deviation re- 
gions, wc do not count the detections of anoma- 
lies occurred in the central deviations region (the 
region enclosed by the solid circle around the pri- 
mary with a radius u s = 0.15). In addition, since 
identifying the existence of both planets and de- 
termining their parameters will be difficult if both 
of the planets' outer deviation regions are located 
at a similar place due to the resulting complexity 
of the interfered pattern of anomalies, we also do 
not count detections if the time interval between 
the anomalies induced by the individual anomalies 
is shorter than O.liE- 

In Figure 2, we present the distribution of 
the absolute probabilities of detecting both plan- 
ets as functions of their projected separations, 
Pabs(bj,bs)- The presented probabilities are the 
values averaged over the random orientation an- 
gles in the range of < a < 2ir. In Figure 3, we 
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Fig. 2. — The distribution of the absolute proba- 
bilities to detect both planets of the example lens 
system with two Jovian-mass planets as functions 
of the projected separations of the individual plan- 
ets from the host star, P a bs{b,j,bs)- The obser- 
vational conditions and the detection criteria are 
described in § 4. 

also present the distribution of conditional prob- 
abilities to successively detect the second planet 
under the condition that the first planet is de- 
tected, P CO nd(bj,bs)- One finds that if the two 
planets are located in the lensing zone, detecting 
both of them is possible with absolute probabilities 
of P a bs 2%. One also finds that once a planet 
is detected, the probabilities to successively detect 
the second one are P CO nd ^ 10%. 

5. Conclusion 

We propose a new channel of detecting extra- 
solar systems composed of having multiple plan- 
ets by using microlensing. In this method, mul- 
tiple planets are detected when the source tra- 
jectory passes both of the outer regions of devi- 
ations induced by the individual planets. From 
the estimation of the probabilities to detect both 
planets of an example Galactic lens system com- 
posed of two Jovian-mass planets around a star 
with ~ 0.3 Mq, we find that if they are located 
within the lensing zone, both planets can be de- 
tected with a non-negligible absolute probabilities 
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Fig. 3. — The distribution of the conditional prob- 
abilities to successively detect the second planet 
under the condition that the first planet is de- 
tected, P a bs(bj,b s ). 

(P a bs si 2%) and a substantial conditional proba- 
bilities (P con d > 10%) of successively detecting the 
second planet under the condition that the first 
planet is detected. The proposed method has an 
important advantage of allowing one to accurately 
determine the planet parameters because the light 
curve anomalies induced by the individual planets 
can be well described by simple single planet lens- 
ing models. 

This work was supported by a grant (2000-015- 
DP0449) from Korea Research Foundation (KRF). 
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